Supercooled confined water and the Mode Coupling crossover temperature 
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We present a Molecular Dynamics study of the single particle dynamics of supercooled water 
confined in a silica pore. Two dynamical regimes are found: close to the hydrophilic substrate 
molecules are below the Mode Coupling crossover temperature, Tc, already at ambient temperature. 
The water closer to the center of the pore (free water) approaches upon supercooling Tc as predicted 
by Mode Coupling Theories. For free water the crossover temperature and crossover exponent 7 are 
extracted from power-law fits to both the diffusion coefficient and the relaxation time of the late a 
region. 
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The effect of supercooling on the dynamics of liquids in 
confined environments is a research field that has become 
more and more popular over the last few years. Out of an 
extremely rich phenomenology that shows diversification 
of specific behavior depending on the size of the particles, 
the confining geometry and the specific interaction with 
the substrate, some general trends can none the less be 
extracted • In fact two competing effects seem to be 
the main contributions to the modification of the dynam- 
ics of the confined liquid with respect to the bulk phase: 
the bare geometric confinement and the interaction with 
the substrate. In particular there is evidence from exper- 
iments that liquid- wall interactions can lead to a layering 
and a decrease of mobility close to the substrate, with a 
substantial increase of the glass transition temperature. 
This effect is stronger for attractive interactions between 
the substrate and the liquid. In some liquids experimen- 
tal evidences show two distinct dynamical regimes Jj]||] . 

Among liquids water plays a most fundamental role on 
earth. The study of the dynamics of water at interfaces 
or confined in nanopores as a function of temperature 
and hydration level is relevant in understanding impor- 
tant effects in systems of interest to biology, chemistry 
and geophysics 

In particular the single particle dynamics of water con- 
fined in nanopores have been studied by different experi- 
mental techniques such as neutron diffraction and nuclear 
magnetic resonance || . A slowing down of the dynamics 
with respect to the bulk phase is observed. Neverthe- 
less the details of the microscopic dynamic behavior of 
confined water upon supercooling are still unclear. 

Below 235 K bulk liquid water enters the so called no 
man's land Q, where nucleation processes, most likely 
triggered by the presence of impurities, take place and 
drive the liquid to the solid crystalline phase, prevent- 
ing the experimental approach to the glass phase ||[|. 



The dynamical behavior of the bulk water simulated 
upon supercooling with the use of the Simple Point 
Charge/Extended (SPC/E) site model potential [To[] , fits 
in the framework of the idealized version of Mode Cou- 
pling Theory (MCT) [jyj, predicting a temperature of 
ideal structural arrest, or crossover temperature Tc, that 
coincides with the so called singular temperature of wa- 
ter p^ ] . This behaviour has been observed several years 
ago by computer simulation |l3] ] and substantiated by 
experimental signatures and further simulation and 
theoretical works [jl5|,^6) . When a liquid approaches the 
crossover temperature Tc MCT predicts that the dynam- 
ics is dominated by the cage effect. After an initial bal- 
listic motion, the particle is trapped in the transient cage 
formed by its nearest neighbours. Once the cage relaxes 
the particle enters the Brownian diffusive regime. Be- 
low Tc, according to the idealized version of MCT, the 
system becomes non-ergodic. In real structural glasses 
hopping processes restore ergodicity. Tc is therefore a 
crossover temperature from a liquid-like to a solid-like 
regime. These activated processes are not relevant above 
Tc for most liquids. 

Until now, no systematic computer simulation stud- 
ies of the microscopic dynamics of confined water upon 
supercooling have been attempted. 

Thermometric studies NMR spec- 

troscopy |]r7|-p0f , neutron diffraction [ l?],|l],^2| and X-ray 
diffraction |23| show evidence that two types of water are 
present in the confining pores, free water which is in the 
middle of the pore and bound water which resides close 
to the surface. Free water is observed to freeze abruptly 
in the cubic ice structure. Bound water freezes gradually 
but it does not make any transition to an ice phase [E3| . 
Layering effects of water close to the substrates have been 
observed in all the simulations for different geometry and 
water-substrate interaction |E|J24| . 
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In this letter we present evidence from Molecular 
Dynamics (MD) simulations that water confined in a 
hydrophilic nanopore exhibits two distinct dynamical 
regimes. In particular the fraction of free water molecules 
behaves similarly to the bulk phase, and its dynamics is 
consistent with several MCT predictions. Power laws fits 
based on MCT yield the crossover temperature of the 
fraction of free water molecules; simultaneously, we show 
that the bound water molecules are already below Tc at 
room temperature. 

In our simulations water has been confined in a sil- 
ica cavity modeled to represent the average properties 
of the pores of Vycor glass j25|. Water- in- Vycor is a 
system of particular interest [ 26[ , since the porous silica 
glass is characterized by a quite narrow pore size distri- 
bution with an average diameter of 40 A. The pore size 
does not depend on the hydration level and the surface 
of the pore is strongly hydrophilic. Moreover the water- 
in- Vycor system can be considered as a prototype rep- 
resenting more complex environments of interfacial wa- 
ter. Different from previous work on the dynamics of 
water close to a planar regular silica surface ]27j we use 
a cylindrical geometry with a corrugated surface. Wc 
have constructed a cubic cell of silica glass by the usual 
quenching procedure. As described in detail in a previ- 
ous work [p5| inside the cube of length L = 71.29 A we 
cut out a cylindrical cavity of diameter 40 A and height 
L by eliminating all the atoms lying within a distance 
R = 20 A from the axis of the cylinder, which is taken 
as the z-axis. After elimination of silicon atoms with less 
than four oxygen neighbours, we saturate the dangling 
bonds of oxygen atoms with hydrogen atoms, in analogy 
to the experimental situation |2(|. 

Water molecules described by the SPC/E model are 
introduced into the cavity. The water sites interact with 
the atoms of the rigid matrix by means of an empirical 
potential model, where different fractional charges are 
assigned to the atomic sites of the silica glass and where 
the oxygen sites of water additionally interact with the 
oxygen atoms of the substrate via Lennard-Jones poten- 
tials 1 25 2^] . The MD calculations are performed with 
periodic boundary conditions along the z-direction and 
the temperature is controlled via coupling to a Berendsen 
thermostat |29|]; the shifted force method is used with a 
cutoff at 9 A to truncate long-range interactions j|o| . 

In this work we consider water at a density cor- 
responding to the experimentally determined level of 
full hydration p6[ . In the chosen geometry this corre- 
sponds to N w — 2600 water molecules and to a density 
p = 0.867 g/cm 3 . We investigate the dynamical behav- 
ior of the confined water for five temperatures, namely 
T = 298, 270, 240, 220 and 210 K. 
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FIG. 1. Mean square displacement (MSD) along 

the z-direction of free water for temperatures (from top) 
T = 298, 270, 240, 220, 210 K at full hydration (N w = 2600); 
in the upper left inset the density profiles at T = 298 K 
(continuous line) and T = 210 K (dot-dashed line) along 
the pore radius are shown; the lower right inset shows the 
diffusion coefficient D vs temperature T for the direction z 
(open circles) and xy (open diamonds). Full lines are power 
law fits to the data, given by D = 9.03(T/185.3 - l) 2 ' 21 and 
D = 7.70(27194.6- l) 1 ' 90 for the z and xy directions, respec- 
tively. D is in cm 2 /s and T is in K. 

In the following we focus on the single particle dynam- 
ics of the water molecules contained in the pore and test 
some of the main predictions of MCT. The radial density 
profile of the water oxygen atoms, normalized to bulk wa- 
ter density, is displayed for T — 298 K and T — 210 K 
in the inset at the upper left corner of Fig. 1. It shows 
the high hydrophilicity of the pore in the form of den- 
sity oscillations close to the substrate. These oscillations 
are not sensitive to supercooling. The average density 
of water molecules at T = 298 KforO<i?<15A 
is p = 0.897 g/cm 3 , for 15 < R < 18 A is p = 1.079 
g/cm 3 and for R > 18 A (a depletion layer) is p = 0.493 
g/cm 3 , where R — y (x 2 + y 2 )- In a preliminary anal- 
ysis done at ambient temperature as a function of hy- 
dration level J3TJ we found that, due to the presence of 
strong inhomogeneities in our system, a fit of the total 
correlators to an analytic shape can be carried out only 
by excluding the subset of molecules in the double layer 
close to the substrate (R > 15 A) which we now identify 
with the so called bound water. In the following we show 
that this shell analysis keeps its validity upon supercool- 
ing. Wc analyze the mean square displacement (MSD) 
and the self intermediate scattering function (ISF) sepa- 
rately for bound water and for the remaining inner layers 
(R > 15 A), which we identify with the free water. Since 
no asymptotic free motion is possible in the xy plane, we 
separately analyze the dynamics within this plane and 
along the pore z-axis. In the main frame of Fig. 1 the 
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MSD of free water is displayed along the non-confined 
z-direction for the investigated temperatures. As super- 
cooling progresses a plateau region due to the cage ef- 
fect develops after the initial ballistic motion, starting 
around t — 0.3 ps. The MSD in the ^-direction (not 
shown) is characterized by a similar trend, but the dy- 
namics is slower. In the low right corner the inset shows 
the behavior of the diffusion coefficient D extracted from 
the fit of the MSD in the Brownian regime region for 
both the z and the xy direction together with fits of 
the MCT predicted power law D ~ (T - Tc) 1 to the 
5 data points. In the z-direction, we obtain Tc — 185.3 
K and 7 ~ 2.21, which are similar to the values found 
for the SPC/E bulk water for ambient pressure, namely 
T c ~ 186.3 K and 7 - 2.29 Hll- In the xy-direction, 
where the dynamics is slower, the critical temperature of 
free water increases to Tc — 194.5 K and 7 ~ 1.90. 
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FIG. 2. Intermediate scattering function (ISF) for free 
water in the zy-direction at the peak of the structure fac- 
tor (Qmax = 2.25 A -1 ) for the five investigated tempera- 
tures. Curves on the top correspond to lower temperatures. 
Full lines are the MD data and long-dashed lines are the fit 
by Eq. ([j]). In the inset the full layer analysis is shown for 
T = 240 K. The central curve is the total ISF, the upper 
curve is the bound water contribution and the lower curve is 
the free water contribution to the total ISF. 

In Fig. 2 we show the ISF of free water at the peak 
of the oxygen-oxygen structure factor along the in- 
direction, Qmax — 2.25 A -1 , as a function of temper- 
ature. The free water molecules inside the pore show, 
similar as in SPC/E bulk water, a diversification of the 
relaxation times as supercooling proceeds. The plateau 
region stretches as Tc is approached. The long time re- 
gion, the so-called late a region, is expected to have a 
stretched exponential decay for a liquid approaching Te- 
la the same figure, the fit of the function 



to the data points is shown, where A(Q) — e~ a Q / 3 is 
the Debye- Waller factor arising from the cage effect with 
a the effective cage radius. r s and 77 are, respectively, 
the short and the long relaxation times, and (3 is the 
Kohlrausch exponent. Obviously the Gaussian form of 
the fast relaxation can be only an approximate one. 

In the inset of Fig. 2 we show the full layer analy- 
sis for T = 240 K as a representative case. The topmost 
curve shows the behavior of bound water, while the lower 
curve is the ISF of free water (identical to the curve in 
the main picture). It is clearly seen that bound water 
is below Tc, since the correlation function does not de- 
cay to zero on the nanosecond time scale. The central 
curve is the total ISF of confined water. It displays a 
strong non-exponential tail, which cannot be fitted by a 
stretched exponential function. Our layer analysis shows 
clearly that the contribution of free water can be sepa- 
rated from the one of bound water and that the stretched 
exponential function is able to give a very good fit to the 
late part of the a region in the free water subsystem as 
we supercool. 
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FIG. 3. The upper part of the figure shows log-log plots 
of the inverse relaxation times, t^ 1 , as a function of (T — Tc) 
along the z-direction (left) and the xy direction (right). The 
full lines are power law fits given by l/r z ~ (T — 185. 3) 2 ' 11 on 
the left and l/r xy ~ (T - 194.6) 1 ' 90 on the right. The lower 
part shows the T dependence of the Kohlrausch exponents j3 
along the z (left) and xy directions (right). 



The values of the Kohlrausch exponents (3 and the re- 
laxation times 77 extracted from the fits of Eq. (0) to the 
ISF data (Fig. 2) are reported in Fig. 3 for free water 
as a function of temperature for the z and the xy direc- 
tions. The T dependence of j3 and 77 is in agreement 
with MCT. MCT also predicts that the inverse of the a- 
relaxation time 77 vanishes with the same power law as 
the diffusion coefficient I/77 ~ (T — Tc) 1 ■ In the upper 
part of Fig. 3 we show log-log plots of the inverse relax- 
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ation times as function of (T — Tc) as points together 
with power law fits as continuous lines. From the fits we 
obtain values very close to the ones obtained above from 
the power law fit of the diffusion coefficients (see inset in 
Fig. 1). This result seems therefore in agreement with 
the MCT prediction that the 7 exponent should be inde- 
pendent of the quantity investigated. In pure SPC/E wa- 
ter the quantity Dt was found relatively constant along 
isochores. None-the-less a slight increase of the product 
was found on cooling most likely due to the progressive 
breakdown of MCT on approaching Tc, see Fig. 5c of 
Ref. p5| . In our case a definitive conclusion on the cou- 
pling of D and r cannot be made at the present stage 
due to the fluctuations in the data. The cage radii ex- 
tracted from the value of the Debye- Waller factor along 
both the xy and z directions range from a = 0.51 A for 
T = 298 K to a = 0.45 A for T = 210 K. The t s values 
are t s ~ 0.2 ps. These values are very similar to those of 
SPC/E bulk water 

In summary, we have presented evidence that the dy- 
namical behavior of SPC/E water confined in a silica 
nanoporc upon supercooling can be analyzed in terms of 
two subsets of water molecules with clearly distinct dy- 
namical regimes, in agreement with signatures found in 
experimental studies on confined liquids |Q-§|. Due to 
the presence of a strongly attractive surface we do not 
find, at variance with MD studies on other confined liq- 
uids, density oscillations all over the confining space Q , 
or a continuous behaviour in going from the center to the 
substrate Q . The bound water being close to a strong 
hydrophilic surface suffers a severe slowing-down already 
at room temperature. Power laws fit based on MCT pre- 
dict that for free water ideal structural arrest would occur 
at Tc 186.3-ftT and 194.5 respectively along z and xy di- 
rections for the full hydration level of the pore. For the 
investigated quantities dynamics appears well accounted 
for by the idealized MCT of supercooled liquids. In this 
respect the predictions of the idealized MCT appear to 
be robust and able to describe also confined molecular 
liquids provided that the effects of the interaction with 
the substrate are properly taken into account. Experi- 
mental evidence of a possible MCT behavior of water-in- 
Vycor have been observed |||3j|. Therefore the analysis 
presented here represents an important step towards the 
understanding of slow structural relaxation of highly non 
trivial glass forming systems. 
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